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A discrete system of coupled waves (with nonanalytic dispersion relation) is derived in the context 
of the spectral transform theory for the Ablowitz Ladik spectral problem (discrete version of the 
Zakharov-Shabat system). This 3- wave evolution problem is a discrete version of the stimulated 
Raman scattering equations, and it is shown to be solvable for arbitrary boundary value of the two 
radiation fields and initial value of the medium state. The spectral transform is constructed on the 
basis of the 9-approach. 
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I. INTRODUCTION 



This paper relates the study of the following discrete coupled system for the three fields Ai{6,n,t), A2{9,n,t) and 
m- q{n,t) 



Ai(e,n,t) - Ai(e,n~l,t) = e"'"''g(n, t)A2(6', n, t) 

A2{e,n,t)- A2{e,n-l,t) = -e"^g(n, t)Ai(0, n, t) (1) 



O' qt{n,t)^^^^ d9e""'iAi^A2){9,n,t) 

• where 9 G [— tt, +7r], n G Z and t G R^. The interaction term here above is defined as the coupling factor 



[Ai*A2){9,n,t) = 

g{e,t)A^{9,n- l,t)A2{9,n,t)+g{9,t)Ai{9,n,t)A2{9,n~l,i) (2) 



where g{9,t) is an arbitrary function in L^([— tt, +7r]) (which could also be time dependent), and where the energy 
I ratio p{n, t) at the site n (it will be shown that this quantity is actually ^-independent) is defined as 



C/2 
> 

, for the following definition of the boundary values Ii{9,t) and l2{9,t) (input data) 



\h{9,t)\^ + \l2{9M 

\A,{9,n,t)\^ + \A2{0,n,tW ^ ' 



/i = lim Ai{9,n,t), h ^ lim A2{9,n,t). (4) 

n — >-4-oo n— * + oo 

One of the main results is the proof that the system (|l|) with data {initial-boundary value problem) 

(?(n,0), h{9,t), l2{9,t) (5) 

is integrable. In particular this work provides the first instance of an integrable nonlinear discrete system with 
nonanalytic dispersion relation. 

An interesting limit of the above equation arises when the arbitrary distribution g{9) goes to a Dirac delta function, 
for instance 5{9). Then the system reads (now Aj{n,t) denotes Aj[9,n,t)\e=Q) 
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Ai{n,t) - Ai{n-l,t) = q{n,t)A2{n,t) 

A2in,t) ^ A2in-l,t) = -q{n,t)Aiin,t) (6) 
qtin,t) = p{n,t)[Ai{n-l,tjA2{n,t) + Ai{n,t)A2{n-l,t)] 

and it is called the sharp line limit of (]l|) . Although a definite physical application of such an equation does not exist 
by now, it still can be understood in the following way. Ai and A2 arc the two envelopes of some high frequency 
(HF) oscillations (say at frequency wi and L02) which interact resonantly on each site n with a medium constituted of 
oscillators of envelope g(n) and frequency ft = uji — 0J2, with a coupling intensity proportional to the relative amount 
p{n) of the total HF energy which has reached the site n. Then the physical data are the input values Ai^2{0, n, t) at 
n = 00 of the HF external excitations, and the initial state of the medium oscillators. 

The method used to build and solve ([^) is the inverse spectral (or scattering) transform (1ST) well known also as 
the nonlinear Fourier transform Indeed the method, in its principle, works like a Fourier transform. It associates 
to the field, solution of a nonlinear evolution equation, its spectral transform which evolves linearly. Then the field at 
time t is reconstructed from the spectral transform at time t by solving the inverse spectral problem. This method has 
been widely studied and extended to various interesting nonlinear evolution problems. We are particularly interested 
in three types of extension which will be used all three together. 

The first extension involved here is the use of discrete spectral problems to solve discrete (in space) nonlinear 
evolution equations. Famous instances of integrable discrete systems are the Toda lattice [|j and a special discrete 
version of the nonlinear Schrodinger equation which has been proposed and integrated by Ablowitz and Ladik by 
using a discrete version of the Zakharov-Shabat spectral problem [|j . This so called Ablowitz-Ladik spectral problem 
and the related nonlinear differential-difference equations with polynomial dispersion relations have been extensively 
studied (see Q - [0). Recently a different version of the discrete Zakharov-Shabat system has been proposed in 
order to keep the canonical Poisson structure of the continuous case Q. 

The second domain considered here concerns the extension of the spectral transform to the case of nonanalytic 
dispersion relations p9| ||20|. T he first instance of such an integrable system is the self-induced transparency (SIT) 
equations of McCall-Hahn ||2l| which was shown to possess a Lax pair in |2|] , was given a N-soliton solution in p^ ] 
and later studied and completely solved in ||2^ . These systems generically describe wave- wave interactions for which 
some boundary value are prescribed. These boundary values are strongly dependent on the physical problem under 
consideration. For instance the problem of superfluorescence in two-level media results in the same equation as SIT 
but with different boundary values and consequently quite different generic properties resulting mainly from a linear 
but non homogeneous evolution of the spectral transform ]2^ ] . 

The third extension used is the generalization of the solution of an evolution equation with a nonanalytic dispersion 
relation to the case of arbitrary boundary values p6t . In this case, the evolution of the spectral transform can be not 
only non homogeneous but also nonlinear and still has interesting physical application. In particular the problem of 
stimulated Raman scattering of a high energy long Laser pulse in a gas has been solved by this technique ||2^] . 

The paper is organized as follows. In sect. 2 we summarize the method of solution of the system (^) and provide 
there only the resulting formulae. 

In sect. 3 the principal Lax operator or, more precisely, the associated spectral problem (a special reduction of the 
Ablowitz-Ladik spectral problem) is used to define the spectral data (or nonlinear Fourier transform). This is done 
by selecting the basic set of Jost solutions and then proving that they obey a Riemann-Hilbert problem in the spectral 
parameter. 

The inverse spectral problem is then solved in sect. 4 by means of the 9-formulation of the spectral problem, 
which means that the potentials are reconstructed from the spectral data. There the compatible reductions are also 
considered, which will allow in particular to obtain simpler integrable equations with an easier interpretation. 

The sect. 5 is devoted to the formulation of the inverse spectral transform on the basis of the 5-problem. More 
precisely, having previously shown that the spectral problem (for the principal Lax operator) leads naturally to a 
i9-problem, we prove here the reverse statement. This is useful in the following for considering the 9-problem itself as 
the starting tool. 

The general discrete integrable systems with nonanalytic dispersion relations is then constructed by requiring a 
time evolution of the spectral transform with a non-analytic dispersion law and a non-homogeneous term. 

These results are used in sect. 6 to prove that indeed the system (|l|) with the arbitrary boundary values (||) is 
solvable. That means that we obtain the time evolution of the spectral data in terms of the boundary values and the 
spectral transform of the initial datum g(n, 0). An interesting case corresponds to the growth of the field q{n,t) on 
an initial medium at rest, that is for g(n, 0) = 0. The method furnishes in such a case the explicit output values of 
the HF fields Ai{9, n, t) and A2{0, n, t) for n — —00. 
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II. SOLUTION OF THE SYSTEM. A SUMMARY 



The general method to generate solutions of ^ is sketched hereafter. The starting point is the spectral transform 
of the initial datum q{n, 0), namely the set of two scalar functions a and /3 defined on the unit circle, a sequence of 
N discrete points kj outside the unit disc to each of which are associated N complex constants Cj . This set is given 
at i = as 

a(C,0), /?(C,0), C = e''; Q(0), kj, |%|>1, j = l,..,iV. (7) 

In the language of the scattering theory, a is called the reflection coefficient, (3 the transmission coefficient, N the 
number of hound states kj and Cj the related normalization coefficients. The effective construction of these data from 
g(n, 0) is displayed in sect. 3, but here we just consider that the set is given and we show how to build from it 
a solution of (|l|). It is worth mentioning that in the linear limit case of small q{x,0), a(C, 0) becomes the Fourier 
transform of q (with parameter 2^), 0) become 1 and all the Cj's vanish (no discrete spectrum, or else no solitons 
in the linear limit). 

The first step is to construct the spectral transform at time t by solving 

d,a^a^{\h{e,t)\' -\h{e,t)\')-{g + g)hh (8) 

1 2 ^,,,2N 



dtkj = 0, (11) 
dtC,{t) ^ -Q(i)^p£^7^(.9 + 5)(|/i|' - 1/2^ (12) 

Unlike in the continuum case, the transmission coefficient j3{C,, t) cannot be computed directly from q;(C, t) and it 
becomes necessary to solve the equation 



(14) 



In the integrals here above, P denotes the Cauchy principal value and C the unit circle in the complex plane. 

Although not elementary, in the case when Ii and I2 are given independently of a and (3 the above system of 
equations is linear and can in principle be explicitly solved as soon as the initial data (|^) and the boundary values 
(^ are known. Before going further, it is already worth remarking that if the quantity /1/2 does not vanish, then the 
evolution for the reflection coefficient a(C, t) has a non-homogeneous term. Consequently the solution can grow on 
the initial vacuum q{n, 0) = which has the spectral transform 

a(C,0)=0, /3(C,0)-1; TV = 0. (15) 

The second step consists in solving the following system of linear integral equation for the unknowns 4>i{k, n, t) for 
1^1 <1 

The solution of (|^) then reads (last step) for C = e'^ 



A2) ' n MO 

q{n+l,t) = -^'f^\n,t) (18) 
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with (j)2 the coefficient of in the Laurent expansion of the solution 02 (fc, t). This achieves the solution of the 
nonlinear system with the arbitrary boundary values (|^) as a sequence of linear operations. 

An interesting information here is the output values (vs. the input values) of the fields Aj (values for n — oo) 
which will be proved to be 



for the input 



A2 I 1 + U|2 I -Iia(3 + 



AA ^ h 

A2 I n^ + co \ I2 



(19) 



(20) 



(note the necessity to compute not only the refiection coefficient a at time t but also the transmission coefficient /3). 
This result, besides having a physical interest, has the nice property of being explicit. Indeed it does not require 
solving the integral equations (|6|) . Actually, when the system (|^) is viewed as describing the interaction of radiation 
components Aj with matter, the relevant (measured) physical information is the output values of the radiation 
components. 

Now we can compute the ratio p{~oo, t) of transmitted photon number, defined in (0), as 



pin,t) ™. i^J^, (21) 



n — > — 00 



while we have obviously 



p{n,t) —> 1. (22) 

n — ^+oc 



This is unlike in the continuous case for which we would find /?(— oo,t) = 1, and results effectively from the discrete 
nature of (|l]) . Indeed a direct calculation leads to the following total photon number non-conservation 

|Ai(C,n- l)p + |A2(C,n- 1)P = (1 + \q{n)\^){\A,{C,n)\' + \A2{C,n)\'). (23) 

As a consequence we obtain from (|^) and the above relation 

00 

p-\n,t)^ n (1 + (24) 

i— n+l 

which proves in particular that the energy ratio p{n, t) indeed does not depend on the variable 9. 



III. THE SPECTRAL PROBLEM 



In the case of the discrete variable, a spectral problem is understood as a difference equation for some unknown /i(ri), 
which involves explicitly an external parameter, the spectral parameter, and a set of given n-dependent coefficients, 
the potentials. Solving a spectral problem results in defining the set of spectral data (functions of the parameter k) in 
such a way that they are in bijection with the set of potentials (in some given class of functions). We shall work here 
in the space of 2 x 2 matrices and adopt an equivalent form of the reduced Ablowitz-Ladik spectral problem used in 
[p[ to integrate the discrete nonlinear Schrodinger equation. In our case we are able to write the spectral transform as 
a i9-problem for the matrix /i(fc, n), which results in a simple formulation of the inverse problem together with a very 
convenient tool for building and solving nonlinear evolutions, in particular those with nonanalytic dispersion relations 
and boundary value data. 

A. Equation and Jost solutions 

Let us consider the discrete spectral problem 

p{k, n+l) - A"^ /i(fc, 7i)A = Q{n+l)p,{k, n+l), (25) 
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with the following definitions 



m 



l/z 




Q{n) 



q{n) 
r{n) 



(26) 



where k is the spectral parameter which belongs to the domain V = C ~ {0, oo}. Up to a change of n —n, 
k 1/k and a rescaling by convenient powers of z of the matrix elements of n it is just the special reduction of the 
Ablowitz-Ladik spectral problem associated to the integrable discrete nonlinear Schrodinger equation. The solution 
of (|25|) possesses the property 



det{^(fc, n — 1)} = det{/i(fc, n)}[l — r{n)q{n)]. 



(27) 



The solution of this spectral problem goes through the construction of some well chosen solutions (the Jost solutions) 
out of some particular asymptotic behaviors. These solutions are denoted by fi^ and are defined by the following 
discrete integral equations 



tiy^{k,n 



Mn(fc," 



Mi2(fc," 
/i22(fc, n 



We will make use also of the notation 



E 

2— n+1 



\ i— — oc 



J 



(28) 



/ 



fc'-"r(i)M+(fc,z) 



n+1 

+00 
■i— n+1 



(29) 



^" '9(*)A*22(fc.* 



.A 



E 

n+1 



r{i)ii^^(k,i) 



(30) 



E fc"- 

i— — 00 

+00 

- E 

t=ri+l 



"*'7(i)M22(^>*) 

r{i)^+^{k,i) 



(31) 



Mil 
M21 



A*2 



± _ / ^^2 



M22 



(32) 



The above integral equations allow then to obtain, for some given class of potentials, the analytical properties of 
the solutions in the domain T) of the complex fc-plane, for all n. The function /i^(fc,n) is holomorphic inside the unit 
circle, the function {k,n) is holomorphic outside the unit circle, the function fj,i{k,n) is meromorphic outside the 
unit circle where it has a finite number N~ of simple poles k~ , the function /zj(fc,n) is meromorphic inside the unit 

circle where it has a finite number of simple poles kj' . Moreover the two solutions /z^ are continuous on the unit 
circle. 



B. Riemann-Hilbert problem and spectral data 

The method to obtain from the integral equations defining /i^ the related Riemann-Hilbert problem is standard. We 
proceed through direct computation of the difference of the two column vectors /x^ and /ij" and obtain the integral 
equation for this difference. It obeys the same equation as the quantity — a~(fc) /i^(fc, n) (the quantity is 
defined below) and, based on the uniqueness of the solution of such equations, we conclude 
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Hi {k, n) — fJ-i {k, n) = —a~{k) fc~" fJ-2 {k, n), |fc| = 1. 
The same approach is appHed to fif and we get 

{k, n) — {k,n) = a'^ (k) fc" 1^.^ {k, n), |fc| — 1 

where the reflection coefficients a^(fc) and a^(fc) are defined (still for |fc| — 1) as 

+00 +00 
a^(fc) = "^('^^ ^^k-'q{i)fi+2{k,i). 

— 00 —00 

For future use we define also 

+00 +00 
r (fc) = 1 -J2q{t)t^2lik,i), = 1 -^r(zK2(fc,*) 



(33) 
(34) 

(35) 
(36) 



which are called the transmission coefficients. Note that, due to the analytical properties of fi2i (resp. IJ.12), fi~ {k) 
can be continued analytically in |fc| > 1 (resp. P~^{k) in |fc| < 1). Actually the vectors /Lt]~ and /x^ have poles where the 
transmission coefficients (i^{k) have poles and we derive from the integral equations (after multiplication hy k — k^ 
and limit k ^ k^) 



Res {k, n)} 
fc- 



(fc-)""Cr + 



Res {^it{k,n)}={ ) (fc+)"C+ 



with the following definitions of the C s 



+00 



+00 



^ q{i)Res {/i2i(fc,«)} 

--\-CO 



i—n-\-l 



/ +00 V 

- E (fc;)"-^9(0Res{^+(fc,z)} ^ 

i—n-\-l 3 

+ CXJ 

- E ''(«)R-es {A*]*2(fc>«)} 

2— n+1 J 



E (fc7)V(z)Re_s{Mri(fc,z)} 



(37) 



(38) 



(39) 



+00 

Ct= E (^.^)"'9(*)Res{^+(fc,z)} (40) 
I— —00 J 

which are called the normalization coefficients 

Since the vectors /i^ and /x^ are holomorphic, we can write down their integral equations evaluated respectively in 
fc^ and kj and compare them to the above integral equations for the residues. We obtain 

Res {Mr(fc,«)} = {kJ^'CJ^x^{kJ,n) (41) 

k~ 
3 

Res {t4{k,n)) - (fc+)"C+A*+(fc+, n). (42) 

Finally the Riemann-Hilbert problem (|3^ ) ( |3^ ) is completed by the behaviors of the solution on the boundaries 
|fc| = and = cx) of P, which read 

^^t{k,n) (J), ^,-(fc,n) ^-^ (j). (43) 
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The vectorial Riemann-Hilbert problem (|3|), (^4|) with singular points given by (^ij), and the boundary behaviors 
(^) constitutes a closed problem which will be solved in the next section. 

The behaviors of ^i{k,n) at large k and of fi2{k,n) at small k will be useful for the following and we define 



Then 



^l^,{k,n) ^ fin). (44) 



fJ.2i{k,n) — > r{n)f{n) (45) 

k — >oo 



and one easily gets that f{n) satisfies the integral equation 



+00 

/(n)-l-^g(z)r(z)/(z), (46) 
n+l 



the solution of which is 

+00 



f{n)^l[[l~q{i}r{i)]. (47) 



ri+l 

The same computation holds for /K^(A;,n) and we have finally 

It can be shown easily by using (|33| ) and ( |3^ ) that the determinant of the matrix ij,{k, n) is analytic in the whole 
domain V. Hence it follows from the Liouville theorem and the boundary values (^), ( p9| ) that 

det{fi{k,n)} = fin). (50) 

Note that, within the reduction r = —q, fin) = p(rt)^^ as given by ([24|). 

C. Asymptotic behaviors and unitarity relation 

By taking the limit at large n directly on the integral equations, the functions obey for = 1 

fi-ik,n) —4 { u-nl^(u\ ?) (51) 



+00 yk-'^a^ik) 1 

M-(^,n) ^ r^?) 1 (52) 



/3-(fc) 



where the following alternative scattering data are defined as 

-\-OQ +00 



a (fc) = '9(*)^22(fc'«)' S+(fc) ^fcV(i)/i+ (fc,i) (55) 

— 00 —00 

+ 00 +00 

;3-(fc) = l-^r(^)Air2(A:,i), ^+(/c) = l-^g(z)Ai2+i(fc,*). (56) 
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The quantities are also called the reflection coefficients to the left (referring to the limit n —oo) when are 
the reflection coefficients to the right. 
It is easy to prove the following relations 



1 — a a+ ' 1 — a a^ 



(57) 



, ^„ /3+= , _ ^ . (58) 
1 — a a+ 1 — a 



Indeed, by using the Riemann-Hilbert problems ( |33| ) and ( p4| ) for (still for |fc| — 1) we have 

— OO — OO 

= a+ (l-^g(*)/iji(fc,i) ) =a+^+ (59) 



+00 4-oc 

/3+ = 1 -^(7(i)At2i(fc,«) +a"^fc"*g(i)^22(fc>*) 

— OO — OO 

^(3-+ a- a-, (60) 

and so on for the other relations. 

Now from (|27| ) the determinant of iJ.{k,n) as n — > — c» can be computed and, by using the behaviors of iJ,{k,n), it 
leads to the relation 

= {l-a-a+)l[[l-r{i)q{i)] (61) 

— OO 

which is called the unitarity relation. 



D. Reduction 



A reduction denotes a simple (possibly algebraic) relation between the potentials (here r{n) and q{n)) for which 
one can derive the counterpart relations for the spectral data. In other word a reduction is a relation which conserves 
the bijection between potentials and spectral data. 

In the case 

r{n) = -q{n) <^ Q{n) = a2Q{n)(72, (62) 

it is easy to check that the function 

i^{k,n) = a2'p{i/k,n)a2, = , (63) 

obeys the same equation as ^{k,n). To compare them it is enough to consider their behaviors as n ^ ±cxd. Since 

a2H (l/fc,n)a2 ^^^^ 1 j ^ ^ 

we conclude that 

a2'p~il/k,n)a2 = fi~^{k,n) (65) 
a-{l/k) ^ -a+{k). (66) 

The same calculation at n — > — oo gives that 

P~{l/k) = P+{k), (67) 
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and also that the alternative scattering data obey similar relations. 
For the discrete spectrum, the relation (p5h implies 



Using the basic relations 



Res {^^t{k,^)} = Res ( ) . (68) 



Res{/(l/fc)} = -(fco)2 Res{/(fc)}, Res{g(fc)} = Res{.9(fc)}, (69) 

fco 1/ko '=0 '^o 



the equation (|68D becomes 



This last equation holds if 



for which it reads 



Res {fit{k,i)} = -{k'^fRes ( 1'^)^''^ ) . (70) 
k+ ^ ' ^ ^ ' \ Aiii(fc,«) ' 



7V+ = 7V-, k-\ = ^ (71) 



Res {/.+ (fc, z)} = -(fc; )^Rcs ( ^^-21 v-^. V ) , (73) 



The above relation together with (|65| ) and ([7l|) implies then 

(73) 

In the case of the reduction (|6^) we will use the following simplified notations (already used in sect. 2) 

a^oT, f3^P', Cj^Cj, kj=k-, N = N+ = N-. (74) 

In order to not over complicate this paper, we do not consider the other reduction r = q for which similar results can 
be easily obtained, but which corresponds to a spectral problem without discrete spectrum. 

Last, it is useful for the following to rewrite the asymptotic boundary behaviors (pT|)-(p4[) within the reduction (and 
with the above notations) and for \(\ — 1 as: 



C "a7 P I n— -00 ^ ' I 1 



M+(C,n) — > ^ " , (76) 



1 + ar 



Similarly, the unitarity relation (|6^) together with the definition ( ^7| ) reads here 

|/3p = (l + |ap)/(-(X3). (78) 
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IV. THE INVERSE SPECTRAL PROBLEM 



The inverse spectral problem consists in reconstructing the potentials q{n) and r{n) from the spectral data 
a±(fc), |fc| = l; 



Cf, k 



\kj\ < 1, J 



l,...iV^ 



\kj\>l, J 



i,..,iv- 



(79) 



A simple way of doing this is to reformulate the analytical properties of the matrix n) in the domain 2? as a 
9-problem. 



A. Inverse problem as a 9-problem 



Indeed, the set of fundamental relations (|33|), (|34|), (41) and (|42|), which contain all the information about the 
analytical properties of fi{k,n), can be summarized in the formula 



d 

^ li{k,n) = fi{k,n)R{k,n), kGD, 
dk 



(80) 



where the spectral transform contains all the information and reads 

R{k,n) = 
/ 



-2i7r 



a+{k)S+ik,l) 
-a-{k)S-{k,l) 

N+ \ 







Y^cpik^kf) 



N- 



y.cj5{k-kj) 



fc-" 
fc" 



The distributions 5^(fc,l) are defined in the appendix and the distribution 5{k) is normalized by requiring that 
JJ dk A dk6{k) — 1. Using the method and tools described in the appendix, and for the behaviors 



y.i{k,n) 



n ) ' f^2{k,n) 



(82) 



(83) 



we have the following solution of this 9-problem 

1 y dXAdX 

Due to the particular structure (^l]) of R{k), the above matrix valued equation has actually to be understood as 
two vectorial integral equations, for /i^(fc,n) and for /i2^(fc,n). As will be seen hereafter, the knowledge of these two 
vectors is sufficient for solving completely the problem. 



B. Reconstruction of the potentials 



The potentials are obtained from the asymptotic expansion of the Jost solutions ^2 a-nd 11'^ of (25) via the formulae 

(84) 



where ^^^2 ^''("■) is the coefficient of fc ^ in the Laurent expansion for fc ^ cxd of Hi2{k, n), and /iJi ^^("') ^fic coefficient 
of k in the Taylor expansion for fc ^ of /i^]^(fc, n). In particular from (p3[) we get 



+(1), 



rin+1) = ^ f dCa-iOC''-^ti22iCn)-J2Ciik-r^-' ^22(^7,^) 
^i" Jc 

1 r 

q{n+l) =2-/ dCa+(C)Cfi+,{(:,n) + Y,C+{k+rfit,{k+,n) 



(85) 
(86) 
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One could easily check that the potentials given by ( Js^ ) and ( p5| ) do obey the reduction (^) when n) obeys (|6^), 
a(fc) obeys (ll), obey (|7|) and obey 0. 

Finally the transmission coefficients f3~{k) and (3^{k) are computed from their definitions (p6|) , where the entries 
lJ.2i{k, n) and ^jLi2{k, n) are obtained from the solution and /i^ by using the explicit relations (pSf), (|3^). Equivalently 
one can use the relations (|5^ ) and the behaviors 

/3+(C) = [l-a+a-] lim ^l22iC,n) (87) 

n— > — oo 

/3-(C) = [1 -a+a"] lim fitiiC,n) (88) 

n^ — oc 

Remark. From the other components of /if and fj.2 we obtain in ( ^5|) the following relations 

oo oo 

f^2t'\n) = E ^'^I'^i + 1)' Mn^'^(") = E + (89) 

n+l n+1 



52) and 



to get (Id = 1) 



V. THE METHOD OF THE 9-PROBLEM 

We have shown in the preceding sections that the spectral problem (|25| ) can be mapped to the i9-problem 
namely 

^li(k)^li{k)R{k), kev, (90) 

with the boundary behaviors (p2|). The solution of such a boundary value problem in the complex plane solves the 
Cauchy-Green integral equation 

The purpose of the following is to show that the above integral equation, for the unknown fi and the datum R, can 
be taken as the starting tool. More precisely we shall show how a parametric dependence of R (on an integer n and 
on a real t) induces the spectral problem (Ea) and a nonlinear evolution equation. 



A. The principal spectral problem 

We restrict this study to off-diagonal matrices R{k) and consider the integral equation (iQll) as the given tool. If 
R{k) depends now on an external integer n, the solution ij,{k,n) solves then the 9-problem po[ ) with the behaviors 

Mfc,n) [l f^\+kf,^'\n) + ... (92) 



^(^-) (ill 



(93) 



where the functions /, /', g and h have to be evaluated. The determinant of ii{k,n) is analytic in V as indeed the 
off-diagonal structure of R(k, n) implies 

|=det{Mfc,«)} = 0, (94) 
ak 

and from the above behavior the Liouvillc theorem implies 

det{M(fc,n)} = /(n) = /'(n). (95) 
We chose now the following explicit dependence of R{k,n) on the discrete variable n 

R{k, n+l) = A(fc)-ii?(fc, n)A(fc), (96) 
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with A(fc) defined in (^6|). The basic fundamental property which allows us to derive from the choice ( |96| ) a difference 
equation for fi in the variable n is the following 

d 

^H{k,n) ^0, i7(fc,7i) = ^(fc,n+l)A(fc)-V(fc,n)"^A(/fc). (97) 
dk 

The above function H{k,n) can then be reconstructed from its behaviors on the boundary of P (A: = c» and k = 0) 
which read from (^) and (93) 

^(f{n)-g{n+l)fig\n) g{n+l) 
fc-o f{n) \ -/(n+l)Ai^VW /("+!). 
H{k,n) ,v'-vpL'H'v ]_ (99) 



fe^oo /(n) V/i(n+l) /(n) - /i(n+l)Ail2^^W 
Since H{k,n) is analytic, these two behaviors are equal, which implies the following four equations 

f{n+l) = ,f{n) - g{n+l)pii\\n) 
f{n+l) = f{7i)-h{7i+l)ti[-^\n) 
g{n+l) = -,f{n+l)y[-^\n) 

/i(n+l) = -/(n+lV^VW, (100) 
which are solved by first defining the potentials as 

q{n+\) = -^x[-^\n), r{n+l) = -^l^^Xn), (101) 

and hence 

g{n) q{n)f{n), h{n) = r{n)f{n), (102) 

with the recursion relation for f{n) 

f{n) = /(n+l)[l - r{n+l)qin+l)] (103) 

of which the solution is indeed given by (|47[). 
Finally, the solution of the (9-problem ( p7| ) reads 

(fc, ,) = \ ( \ 9(^1) ) (104) 

^ ^ 1 - r(n+l)g(n+l) 1 / 

which can be written with ( p7| ) as the discrete spectral problem 

^i{k, n+l)A{k)-^fi{k, n)-^A{k) = [1 - 0(n+l)]^i (105) 

where 



r(n) 

It is convenient for the following to define the quantity 

U{n)^[l-Q{n+1)]~^ (107) 

such that the equation for /i(fc, n) reads 

f^{k, n+1) = U{n)A{k)-'^fi{k, n)A{k) (108) 
which is precisely the spectral problem (Eq). 
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B. Nonanalytic Dispersion Relations. A Theorem 



We consider now that R{k, n) depends also on an external real t and address the problem of computing the expression 
of the time dependence of the solution /i of (|9|) . The result can be stated as a theorem. 
Theorem. When the spectral transform R(k,n,t) evolves according to 



Rt{k, n, t) = [R{k, n, t) , n{k, t)] + M{k, n, t) 

where 

M{k, n+l,t) A{ky^M{k, n, t)A{k), [A(fc) , n{k, t)] = 0, 
and where fl{k,t) is the nonanalytic dispersion relation 



ZlTT 

the potential Q obeys the following evolution 

Qt{n+l,t) = 

where 



dXAdX dn{X,t) /fc/A 
A -A: a A V 1 



C^3 



2in 



d\ A dX 
2A 



T{X,n,t) 



T{k, n, t) = A{k)-^fi{k, n, t){M{k, n, t) - ^^^^£^}A(fc) n-\k, n+1, t). 

dk 

This theorem is proved hereafter. 
The auxiliary spectral problem 
Let us define the matrix 

n, t) = {lJ.t{k, n, t) — ^{k, n, t)n{k, t)}fi^^{k, n, t), 
and compute its 9-derivative which, from (|8^), ( |109| ) and ( |llC| ), obeys 

^Ki^ ^ ^(fc, n, t){M(fc, n, t) - MM) t). 
To solve the above 9-problem we need the behaviors of Vi as fc — > and V2 as /c — > 00. Since 



n, t) 



^{k, n, t) 



k^o V 



k—>oo 



(109) 
(110) 

(111) 

(112) 
(113) 



(114) 
(115) 

(116) 
(117) 



it is easy to obtain, thanks also to the choice (111) (it would not be so in the case of a regular dispersion relation) 



Viik,n,t) 



V2ik,n,t) 



Consequently the solution reads 



V{k,n,t) = — JJ^^—^ S{X,n,t) 



(k/X 
I 1 



where we have defined 



S'(fc, n, t) — fi{k, n, t){M{k, n, t) ^{k, n, t). 

dk 



(118) 



(119) 



(120) 



Now, with the above value of V, the definition (114) can be written as the auxiliary spectral problem 
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fXt{k, n, t) = V{k, n, t)fi{k, n, t) + n{k, n, t)rt{k, t), (121) 
Fundamental property of V{k, n, t) 

For simplicity of notations, wc omit from now on the variable t. By direct computation the matrix S obeys the 
equation 

S(k, n+l)U{n) = U{n)A{k)-^S{k, n)A(fc). (122) 

The next step consists in seeking an analogous property for V{k,n), by computing the quantities U{n)^^V{k,n+l) 
on one side and A{k)~-^V{k,n)A{k)U{n)~^ on the other side. Using 

A{X)-^S{X,n)A{X)U{n)-^ (^''(^^ ^ A{k)-^ S{X,n) (^''(^'^ A{k)U{n)-^ 

^ A_fc / -g(n+l)sii(A,n) + Asi2(A,n) 

A \r{n+l)s22{X,n) - {l/X)s2i{X,n) 

we obtain finally the required property of V{k,n) 

U{n)-^V{k, 71+1) - A{k)-^V{k, n)A{k)U{n)-^ = P{n) (123) 

T(fc, n) = A{ky^S{k, n)A{k)U{n)-^ = U{n)-^S{k, n+1). (125) 



We have used here above (|122| ) to rewrite T in a more convenient form, and finally from the definition ( |120| ) of S{k, n) 
and the spectral problem ( ^51) it reads 

T{k, n) = A(fc)" V(fc, n){M{k, n) - ^^^}A(fc) fi-\k, n+1). (126) 

dk 

Note that the matrix T(fc, n) obeys a property similar to S{k, n) since it can be checked directly that 

U{n - l)r(fc, n-1) = A(fc)T(fc, n)U{n)A{ky\ (127) 

The evolution equation 



Th e no nlinear evolution of Q{n,t) is now obtained in the usual way by requiring the compatibility between ( IPS 



and (121) which reads 

d . . d 



^A^(fc,n+l,t) = ^{C/(n)A(fc)-V(fc,^)A(fc)} 

= V{k, n+1, t)^{k, n+1, t) + ^i{k, n+1, t)n{k, t). (128) 



By means of (123) it is then easy to obtain the equation 

Ut{n) = U{n)P{n)U{n), (129) 
which readily gives the evolution ( |112| ) since U{n) = [1 — Q{n+1)]^^ . This ends the proof of the theorem. 

C. Reduction 

If we consider the reduction r(n) = —q{n) the matrices 

0(M)^ (130) 



must be compatible with the preservation of the structure of R{n) in the time evolution equation (109). We choose, 
therefore, lui and tU2 analytic inside and outside the unite circle with limit values on the two sides of the circle satisfying 
the symmetry properties (C — e'^) 
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c^r(C,t)=u;2+(C,t) (132) 

W2"(C,t)=^(C,t) 

and we choose mi and m2 as 

m^{k,t) ^m-{k,t)5^{k,l) = to(C, <)(5" (fc, 1) (133) 
m^(k,t) =m+(fc,i)(5+(fc,l) EETO(C,t)(5+(A:,l) 

where 

^ \k\ 

and m((^,t) is a given function defined on the circle |fc| = 1. Note that the discontinuity of uji 

piC,t)=uj+{C,t)-Lu^{C,t) (134) 
is related to the discontinuity of W2 by the formula 

io+{C,t)-co^{C,t) = -piC,t). (135) 
The analytic properties of cui and uj2 are summarized by the formulae (C = k/\k\) 

^=p(C)5(fc,l) (136) 

^ = -p(C)5(fc,l) (137) 

where the distribution 6{k, 1) is defined in the appendix. Requiring that loi —^ for k and UJ2 for k — > cx), 
ri(fc) is defined by the following Cauchy-Green formula 

It results that 

uJi{k)=LJ^{l/k) (139) 



in agreement with the conditions on the boundaries ( 132 ). 

It can be shown that with the choices indicated in (132) and (133) the evolution equation ( |112| ) is compatible with 
the reduction r = —q. This evolution reads 



1 



where 



+^2 (C)CMr2 (C, "-)Mri {Qin+l)-uj+ (C)CM]*2 (C, (C, n+l) 

+a{C) (c^r (0 - ^2 (C))C->r2 (C, n)^l-^ (C, n+l) 
+a(C) {o^t (C) - (C))C"+ V^i (C, n)^Jitl iCn+l) 

-m(c)r'vr2(c,n)Mr2(c,'^+i) + m(c)c"+V^i(c,«)A*^i(c,«+i). (wi) 

with C = e*^ (remember that a = a" and a~''(C) = — a(C)). Note that in computing TiX^n) in ( |112| ) the term 
containing ri(/c, i) must be written as follows 

A-V(ri)-=AA*-i(n+l) = 
afc 

-1 (A- V('^)^^AM"^ ("+!)) - A-V(?^) [-R("), ^] I^H'\n+\) (142) 

which is a well defined local formulation of a 9-problem for a sectionally holomorphic function. 
It is convenient, by using equations (|3^ ) and (^4|), to rewrite ( |l4l| ) in terms of /ij^ and /^^ 

7(^^, n) = p(C)Mii(C, n)A^r2(C, «+l) + |5(C)CMr2(C, n)^lUQ, n+l) 

-m(c)r'vr2(Cn)Mr2(c,«+i) + m(c)c"+V^i(c,n)/^^i(c,«+i) (143) 

which shows explicitly that the evolution equation depends only on the discontinuity of ^l[k) on the unit circle. 
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VI. INTEGRABLE DISCRETE INITIAL-BOUNDARY VALUE PROBLEM 



By using the tools previously developed we prove now that the nonlinear system (|l|) is integrable when it is related 
to the initial boundary value (||). We rewrite hereafter this system in the variable C, — e*^ and with the relation ( |2^ ) 
as 



{n,t) n (l + l9(*,<)n = ^ / 

i— n+1 ^ 



d9e 



\A,^A2)iC,n,t) 



Ai (C, n, t) - Ai (C, n-l,t) = C'"q{n, 0^2 (C, n, t) 
A2 (C, n, t) - A2(C, n-1, i) = -Clin, t)Ai (C, n, i) 



(144) 

(145) 
(146) 



where the interaction term is defined as 

(Ai * A2)(C, n, t) = g{e, t)Ai(C, n-1, i)A2(C, n, t) + 5(0, OAi(C, n, i)A2(C, n-1, i). (147) 
Theorem. With the datum of the initial value q{n, 0) and the following arbitrary boundary values as n +00 

Ai(C,n,i)->/i(C,i), A2{C,n,t) ^ hiC^t). (148) 
the above system is solvable by the spectral transform method. 

A. Proof of integrability 



The proof is performed by showing that the evolution ( |144| ) can actually be written under th e for m (|140|). This 
then gives a unique definition o f the functions m{k) and uj{k) +u;{k) for which the two equ ation s ( 140) and ( |144| ) are 
identical. Hence the evolution (109) of the spectral transform R{k,t) is uniquely given via (132) and (133). 

The first useful property is that the following 5 vectors 



C-"^2(C,^) 



t^tiiCn) 

Ai2l(C,"-) 



M22(C,") 



solve the equation (146). Then, by comparison of their asymptotic behaviors as n 
we get 



(149) 

-00 given in (|l48|) and in 



r"^2(c,«) 



= h 



«)(:ll;:!)+«<K-"(:i:;:! 



Next, to compute the product {Ai * A2){C,n) we make use of the Riemann-Hilbert relations 

/i2-(C,«) - A*+(C,n) = -a+(C)C" ^^t{C,n) 
and rewrite it in terms only of /i^ and 

C"+i(Ai * A2)(C,n+l) = - [g\h\'-g\l2\^] /i+(n)^r2(«+l) 

-[5|/i|'-5|/2|'] CMr2("Ki(«+l) 
-ig + g)hl2C''tJ-l2Hf^i2{n+l) 
+ (g + 5)W"+VfiWA^^i('^+l) 



(150) 



(151) 
(152) 



Then, thanks to the expression Hm), the two equations (|140|) and (144) are identical if and only if 



p(C,i) = -g{0,t)\hi0,t)\' +gie,t)\l2i0,t)\^ 

rnict) = {g{e,t) + g{e,t))h{e,t)i2ie,t). 



(153) 



(154) 
(155) 



Finally the theorem is proved and it remains to compute the evolution of the spectral transform. 
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B. Evolution of the spectral transform 



Evolution of a(C, t) 



The time evolution of R{k, n, t) is given by ( |109| ) with M and 51 defined in ( |133| ) and in (138). Taking into account 
the structure (pl|) of i?(fc, n, t) we have 



ata(C, i) = [c^r (C, i) - ^2 (C, i)] a(C, " (C, t) 



where from ( |138| ) and the Sokhotski-Plemelj formula wc have 

wr(c,i)-c^2"(c,t)^ 



1 

1 

1 



cC'-C 

dC 

cC-C 
dC 



2711 Jq (' — kj 

dC 



2iTi Jc ~ k. 



and the functions p(C,i) and •m{C,,t) are given in (154) and (155) 
As a result the evolution equation of a and Cj can be written 



{g + g)hl2 



/cC'-C 



{9 + 9){\h?-\h?) 



dtC,{t) 



Evolution of P{C,,t) 

The definition (^8^^ allows to obtain readily 



(k 



(156) 
(157) 



(158) 
(159) 
(160) 
(161) 



(162) 
(163) 
(164) 
(165) 
(166) 



l + |a|2' 

which actually can be understood as a Riemann-Hilbert problem on the unit circle. Its solution reads 



(167) 



13+ dt\l + |a|2 



l+\a\ 



1 



27ri JcC,' -kl + |a(C')l" 



(168) 



|fc|<l 



/3 



dC k {\a{C)?)t 



27ri JcC ~k C 1 + |a(C')P' 
Hence, writing the above equation for k = C(l — 0), we get the evolution (jl^). 



(169) 
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C. Evolution of the spectral transform from the Lax pair 



For completeness, we rederive hereafter the preceding formula (evolution of a and in the abse nce o f bound states 
for simplicity), by using the traditional approach for which the starting tool is the Lax pair ( |2^ ) ( 121 ). The met hod 
consists simply in evaluating the asymptotic boundary values as n — > ±cxd on the auxiliary spectral problem (121), in 
which (forget for a while the (n, t)-dependence) 



V{k) = 

ZlTT 



^-^MA)(^Af(A)-^j/. (A)(^ Q ^ 



By using the identity (142), the equation (121) can be more conveniently written as 

/it(fc, n, t) — X{k, n, t)ii{k, n, t), 

where 



ZlTT ,1,1X1 



^^^(A, n, t) {M{\) + [R{X), n{X)]) ^^-\X) J 



(170) 



(171) 



(172) 



By inserting in the above equation the explicit form s of R{ k, n , t) given in (|l|) with (H) and no boimd states, of 
fl{k,t) given in (13J), and of M{k,n,t) given in (131) with (133), we get finally 



(173) 



with the following definition 



X(C) = Cirn + {c^t ^})^] ( ) 

A-nr™ (, ,- , ( A^i2/"22 ^(^12) 

\ IM22J M12P2 



2 
^^22 



(174) 



The main tool is now the asymptotic boundary behaviors ( [75[ ) and ( |76| ) of /i^ which allows to obtain, by taking 
the limit as n ^ +c>o of C"c)(/i^]^(C, n, i), the relations 



at{k) 



[to — [ijJi — )a] (fc) + 



n+l 



= 



1 f dC fk 
lim P (B — - — 

n^oo 27ri / C ^ ^' VC 

-i [TO+(a;+-^2+)a] (/c) + 



[to - (wi - ^2 )a] (C) 



1 



lim - — :P ^ 7 1 , 



^^^"^TO+(C.+ -C.2+)S](C). 



Consequently, with the formula (see appendix) 

1 



lim -^P 

n— »cx3 2ni 



\k\ 



(175) 



the preceding relations result precisely in the required evolution ( |156| ). 

In the same way, by taking the limit as n — > —00 of 9t/ij~]^(C, n, t), we obtain readily the required evolution (Q) of 
the transmission coefficient /3(C, t) (note that there, one should use also the unitarity relation 



D. Time evolution of f{n) 
It could be useful to have also explicitly the time evolution of the quantity /(n).From (BOh we have 
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Mn) = f{n)tT{^it{k,n)fi-\k,n)} (176) 
and then using the auxihary spectral problem ( |171| ) 

ft{n) = fin)tr{X{k,n)}. (177) 
From the expression (|l73| ) we obtain that the trace of X{k, n) is fc-independent and reads 

tr{X(fc,n)} = Y ■f^""^^ii^2i - C^mfi^^fi^^ + (p + p)fJ- 12 fJ-ti 

-{UJ+ - tJ+)/^+ ^+ + (tjf - t^2")/ir2Ai2"l} • (178) 

Due to the analyticity of the function 

^(Wl(fc) - L02{k))Hl2{k)^l2l{k) 

inside and outside of the circle the last two terms in the r.h.s. vanish and we obtain for the evolution equation of /(n) 

p+TT 

ft{n)^ d9 {Cm4l^J'2l-C''mfl^2^J■22 + iP + P)^J'l2fJ'2l} ■ (179) 



This result can also be expressed in terms of the physical quantities Aj and Ij by inverting (150) to get on the unit 
circle 



f^ii ^ M22 



Aili + A2I2 



|/lP+|/ll2 



M2I- A^12- C |J^|2 + |J^|2 

and then by inserting these formulae and those for p + p and m in (|154[) and (|155|) into (|179|). We obtain finally 



/.W^ir..(,, + B™^™. (ISO, 



2 



Note that f{n) is conserved if g is pure imaginary. 



APPENDIX A: MATHEMATICAL TOOLS 
1. Basic distributions 

The distributions i5^(A, 1) and 5(A, 1) have support on the circle C of radius 1 in the complex A-plane and are 
defined by the following formulae 

'''' d\ A dA (5±(A, 1)/(A) = / dC / ((1 T 0)0 (Al) 

■D Jc 

dXAdXS{\,l)f{\)= (f dCfiC) (A2) 
V Jc 



or, equivalently, by 



d\AdX6^{X,l)f{X)^i d9e''^f {{lTO)e''') (A3) 

J — 7T 

dXAdX6{X,l)fW = i f d9e'^ f {e''^) . (A4) 



Note that the distributions ^^(A, 1) can operate on functions which have defined left or right limit on C, while the 
distribution S{X, 1) can operate on functions continuous on C or with support on C. 
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By complex conjugation of the equation (A3) and by the change of variable 6 —6, we obtain 



// dAA(iA<5±(A,l)/(A) =^ / dOe'" f {{1 T 0)e-''^) . (A5) 

JJV J -TV 

Next, by means of the change of variable A — > 1/A (remember that the domain T) docs not contain the point A = 0) 
and — > —9, we obtain 

JJ^d\AdXS^ij,l)f{\)=i d9e''f{il±0)e-^'), 

and consequently 



(A6) 



^A 

Now, through the change of variable A — > A, we obtain 

// dXAd\d^(X,l)f{X)=i f dOe'" f {{1 T 0)e'' 

JJV J-TV 

which implies the second symmetry property 



<5±(i,l) = jT(A,l). (A7) 



(A8) 



5±(A,1) = (S±(A,1). (A9) 

These two relation naturally leads to 

,5+(A,l) = <5-(i,l). (AlO) 
A 

Similar symmetry properties can be obtained for S{X, 1). 

2. The generalized 9- formula 

Let F+ e C'^(D+) where D+ is the open disk of radius 1 centered in the origin and F G C^(CZ?+), let and 
satisfy the Holder condition on the circle C of radius 1 and let F~ vanish at large z. Then by noting F{z) the 
function defined as F~^(z) for z € £)+ and as F^(z) for z € D~ = CD+ we have 

F{z) = 7r-f ^7 —'^^ + ^ -^—dXAdX + 

1 [[ dF/dX^^^- ^^^^^ 



2iTi J Jd- X — z 

where the circle C is anti clockwise oriented. 

If we define the 9-derivative of a function F{z) discontinuous on C as follows 

8F 8F BF z 

— ^{F+{0-F-{0)6{z^) + 4>D.{z)^ + <t>D-{z)-^, (A12) 



where 4>a{z) = 1 for z E A and (jiAiz) = for z ^ A the generalized 9-formula (All) can be rewritten as 



F[z) = !-.((_ ^-p^dXAdX. (A13) 



Stti ././"5+uD- ^-z 



Formula (A12) can be considered as the local formulation of the generalized Cauchy-Green formula (All) 
Subtracting formula ( |A11 ) at z = a we obtain 
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1 /■/■ dF/dXfz-a\ ,^ - 1 /■/■ dF/dXfz-a\ ,^ - 

" ' ' ' dA A dA + — / / — ^ UA A dA (A14) 



27ri y Vd+p a — z \A— a y 27ri J J^,- A — z \A— a 

if a G and an analogous formula if a G . The second integral on the right and side is obtained first by 
computing it on the set Da.e = {A : A G Z?, |A — a| > e} and then by taking the limit e ^ 0. Note that the formula 
remains valid also if F{z) is going to a const ant d ifferent from for z oo. 
If for z ^ oo F{z) F{oo) we can apply ( All ) to F{z) — F{oo) getting 

Fiz) = Fioo)^l-.[Em-I^dC^ 



27ri Jq C ^ ^ 

1 ff dF/dX^^ - 1 ff dF/dX^^ - 

' dXAdX+— dXAdX. (A15) 



27ri J J JJ+ X — z 2ni J J j^- A 

Finally let us note that the Sokhotski-Plemelj formula on the circle reads 



/(C') = ±*^/(C) + ^/7^^ /(C), ICI = 1, (A16) 



/cC'-(itO)c'^" ' JcC~C 

where denote the Cauchy principal value integral. 



3. Limits of integrals 

Let us proof that 

Under the following successive changes of variables 

(^ = 6"', k = e''^, a = d~if, x = na (A18) 

we derive 

1 f df /k\^ 1 ri-^-v)" -ix , , , 

i^/A 7 *(C) = i^^/^ , (A19) 
and taking the limit, for — tt < (p < tt, 

■ -2m J C~k\CJ 2711 J^^ X 
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